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LONGITUDINAL AND LATERAL SPECTRA OF TURBULENCE 
IN THE ATMOSPHERIC BOUNDARY LAYER 

SUMMARY 


In many types of space vehicle and aircraft response calculations, 
engineers employ the Fourier transform to solve the equations of motion which 
describe and predict the ultimate response of space vehicles to the natural 
environment. Thus, the input functions which describe the turbulent character 
of wind loads must be specified in terms of spectra. To develop an engineering 
boundary layer model of the longitudinal and horizontal lateral spectra of 
turbulence, a 150-meter meteorological tower was erected at Cape Kennedy, 
Florida. The tower is instrumented at the 18-, 30- , 60- , 90- , 120- , and 
150-meter levels with wind speed and direction sensors which measure the 
horizontal components of the wind. Temperature sensors are located at the 
18-, 30-, 60-, 120-, and 150-meter levels. Fifty one-hour cases of turbulence 
that occurred during a variety of meteorological conditions were selected for 
analysis. 

Longitudinal (u) and horizontal lateral (v) velocity fluctuation 
spectra were calculated for each level, resulting in a net total of 300 spectra 
for each velocity component. It was found that ^nS(n)/u^, f ] - coordinates, 

or rather Monin coordinates, fail to produce a vertical collapse of the spectra, 
where S (n) is the longitudinal or lateral spectral energy density at frequency 
n (Hz) , u^ ^ is the surface friction velocity (calculated with mean wind and 

temperature profile data) , and f = nz/ u, u is the mean wind speed at height 
z. However, it appears that Monin coordinates will collapse spectra with 
various turbulence intensities at any particular level in the vertical. 

To produce a vertical coUapse of the spectra, it is assumed, for 
engineering purposes, that the spectra in Monin coordinates are shape- 
invariant in the vertical. This hypothesis, which seems to be reasonable, 
permits a practical approach to developing an engineering spectral model of 
turbulence. The longitudinal and lateral spectra for the unstable and neutral 
boundary layers can be represented by the form 



Cf 


nS(n) _ 


m. 


1+1.5 


(f-) 

m 


1 5/3r 


where C and r are functions of stability. In the neutral boundary layer, 

and are proportional to and and in the unstable boundary 

layer they are proportional to and The dimensionless wave 

number associated with th'e peaks of the logarithmic spectra is an 

increasing function of z. In neutral air, f and f are proportional to 
® mu mv 

z and z®*®®, while in unstable air, they are porportional to z®*®^ and z®*^^. 

The inverse Fourier transforms (correlation functions) of these spectra 
are calculated numerically. TOe integral scale L» Is proportional to 

SO that, in the neutral air, L* is a constant and L^' is proportional to z®*^^, 

and in unstable air they are proportional to z®*^® and z®*®®. 

An analysis of the logarithmic spectrum in the inertial subrange, at the 
18-m level, implies that the local mechanical and buoyant production rates of 
turbulent kinetic energy are balanced by the local dissipation and energy flux 
divergence, respectively. In addition, the spectral model implies that the 
dimensionless dissipation rate kj z e/ u®^ ^ is proportional to z®*®®® and z®*®® 

in the neutral and unstable boundary layers, where ^ and e denote 
von Karman’s constant and the dissipation rate. 


INTRODUCTION 


To determine the response of space vehicles, aircraft, tall structures, 
etc. , to atmospheric turbulence, the engineer requires specific information 
about the spectral nature of atmospheric turbulence. This results because the 
equations of motion of these vehicles or structures are linear and are solved 
with Fourier transform techniques and thus the environmental forcing functions 
must be represented in terms of spectra. Motivated by this requirement, we 
have developed a model for the longitudinal and horizontal lateral spectra of 
turbulence for the Kennedy Space Center (KSC) . The longitudinal and lateral 
components of turbulence are the wind fluctuations parallel and normal to the 
mean wind vector (Fig. 1). 
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Figure 1. The relationship between the quasi-steady and the 
instantaneous wind vectors and the longitudinal and lateral 
components of turbulence. 


THE NASA 150-m METEOROLOGICAL TOWER 


To obtain micrometeorological data representative of the Cape Kennedy 
area, especially in the vicinity of the Apollo/Saturn V launch pads, a 150-m 
meteorological tower was constructed on Merritt Island at KSC. The tower 
facility, discussed in detail in a report by Kaufman and Keene [1] , is only 
briefly described here. 


Terrain Features 

Figure 2 shows the location of the facility with rejspect to the Saturn V 
space vehicle launch complex 39. The tower, located about five kilometers 
from the Atlantic Ocean, is situated in a well-exposed area free of near-by 
structures which could interfere with the air flow. 
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Figure 2. NASA Launch Complex 39, Kennedy Space Center, Florida. 


The aerial photograph (Fig. 3) of the terrain surrounding the tower 
(point T) was taken at 1066 m above mean sea level. In the quadrant from 
approximately 300 deg north azimuth with respect to the tower, clockwise 
around to 90 deg, the terrain is homogeneous, being covered with vegetation 
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Figure 3. Aerial view of the terrain surrounding the NASA 150-m 

meteorological tower. 


about 0.5 to 1.5 m high. Another homogeneous area with the same type of 
vegetation occurs in the 135- to 160-deg quadrant. The areas A (230-300 deg) , 
B (90-135 deg) , and C ( 160-180 deg) are covered with trees from about 10 to 
15 m tall. The distance from the tower to areas A or C is about 200 m, and the 
distance to area B is about 450 m and the vegetation ranges from one-half to 
one and one-half meters, as in the area to the north of the tower. To the 
south-southwest, in the 180- to 230-deg quadrant 225 m from the tower, there 
is a body of water called Happy Creek. 



Instrumentation 
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The complete tower facility comprises two towers, one 18 m. and the 
other 150 m high (Fig. 4) . The levels on both towers are instrumented with 

Climet (Model Cl-14) wind sensors. 
Temperature sensors, Climet (Model 
Cl-016) aspirated thermocouples, are 
located at the 3- and 18-m levels on 
the small tower and at the 30-, 60- , 
120- , and 150-m levels on the large 
tower. Foxboro (Model F-2711AG) 
dewpoint temperature sensors are 
located at the 60- and 150-m levels 
on the large tower and at the 3-m 
level on the 18-m tower. Wind speed 
and direction data can be recorded on 
both paper strip charts and analog 
magnetic tapes with an Ampex 
FR-1200 fourteen-channel magnetic 
tape recorder which uses a 14-inch 
reel. The temperature and dewpoint 
data are recorded on paper strip charts. 
To avoid tower interference of the flow, 
the large tower is instrumented with 
two banks of wind sensors. The details 
of how and when one switches from 
one bank of instrumentation to the other 
is discussed by Kaufman and Keene [1] . 
During a test in which the wind data are 
stored on magnetic tape, only one bank 
of instrumentation is used, thus avoiding interruption of the wind data signals 
within any magnetic tape recording period. This continuity of signals prevents 
data processing difficulties when converting analog tapes to digital tapes. 
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Figure 4. Location of instrumen- 
tation on the NASA 150 m meteoro- 
logical tower at Kennedy Space 
Center, Florida. 


Surface Roughness Length (Zq) 

In an earlier report, Fichtl [2] discussed the surface roughness length 
configuration associated with the NASA meteorological tower. This analysis 
was based upon wind profile laws that are consistent with the Monin-Obukhov 
similarity hypothesis. The calculations of were based on wind data 
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obtained at the 18- and 30-m levels and on temperature data obtained at the 
18- and 60-m levels. Most of the measurements were obtained during the 
hours of 0700 and 1600 EST, and the gradient Richardson numbers at 23 m 
(geometric mean height between 18 and 30 m) for the 39 cases ranged between 
-5. 82 and +0. 079. The results of these calculations, shown in Figure 5, show 
the effect the terrain features have upon the surface roughness. Later, an 
analysis of the energy budget at 18 m will show that these roughness lengths 
are too large. 


Zq (meters) 



Figure 5. Tentative azimuthal distribution of the surface roughness 
length at the NASA 150-m meteorological tower site. 


COMPUTATIONS AND INITIAL SCALING 

To establish a spectral model of turbulence for the Kennedy Space 
Center, approximately 50 cases of turbulence were analyzed. The procedure 
used to calculate the longitudinal and lateral components of turbulence 
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consisted of (1) converting the digitized wind speeds and directions (10 data 
points per second) into the associated north-south and east-west components 
and averaging these components over the duration time of each test, (2) cal- 
culating the mean wind speed and direction with the averaged components, (3) 
projecting the original digitized data onto the mean wind vector and subtracting 
the mean wind speed to yield the longitudinal components of turbulence, and 
(4) projecting the original digitized data onto a normal-to-the-mean wind 
vector to obtain the lateral components of turbulence. Trends contained within 
the data were removed by fitting the longitudinal and lateral components of 
turbulence to second order polynomials and, in turn, subtracting these 
polynomials from the component time histories. To reduce computation time, 
the data, with trend removed, were block-averaged over half-second intervals. 
The longitudinal and lateral spectra were calculated by using the standard 
correlation Fourier transform methods given by Blackman and Tukey [3J. 
These spectra were corrected for the half-second block-averaging operation 
with the procedure given by Pasquill [4] and for the response properties of 
the instrumentation. 

To combine the spectra for each level on the tower, it was assumed that 
the similarity theory of Monin [5] for the vertical velocity spectrum could be 
applied to the longitudinal and lateral spectra, so that 


= F(/,Ri) , (1) 


where nS(n) is the logarithmic longitudinal or lateral spectrum associated 
with frequency n(Hz), and u, ^ is the surface friction velocity, or rather, 

the square root of the tangential eddy stress per unit mass. F is tentatively 
a universal function of the dimensionless wave number f and the gradient 
Richardson number Ri. The dimensionless wave number is given by 


nS(n) 


u“ 


o 


nz 

u(z) 


( 2 ) 


Because the tower did not have the capability to measure vertical 
velocity fluctuations, the Reynolds stress, and hence u^ , cannot be calculated 

with first principles; viz, u^^ = -u*w', where u’ and w* are the longitudinal 
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and vertical velocity fluctuations and the overbar denotes a time-averaging 
operator. However, an estimate of the surface friction velocity can be 
calculated from mean wind and temperature profile data. 

According to Lumley and PanofslQ^ [6] , the mean wind profile in 
approximately the first 30 m of the atmosphere is given by 


u(z) 



(3) 


where is von Karman's constant with numerical value approximately 
equal to 0.4, and is a universal function of z/L'. L* is a stability length 
given by 


L' 


u 


»;<o 


du 

dz 


T 



(4) 


where T and 6 are the Kelvin and potential temperatures associated with 
the mean flow. The quantity z/L* is related to the gradient Richardson 
number 


Ri = 


M. ^ 

6 dz 



(5) 


through the relationships 


z _ Ri 

L* ~ (1 - 18 Ri)V^ 


(Ri < -0. 01) 


( 6 ) 


z 

L* 


= Ri 


(-0.01 < Ri ^ 0.01) 


(7) 
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and 


z _ Ri 
V “ 1 - 7 Ri 


(0,1 s Ri> 0.01) 


( 8 ) 


Equation (6) is a form of the KEYPS [6] equation. The function ’t(z/L') 
associated with equations (7) and (8) are given by 


fr 




(-0.01 < Ri:< 0.01) 


(9) 


and 






(0.1 ^ Ri> 0.01) 


( 10 ) 


Lu33oley and Panofslqr [6] have graphically indicated the function 'i' (z/L*) for 
Ri < -0. 01 and the function 



0.044 



1.0674 - 0.678 In 



(Ri < -0.01) 


( 11 ) 


faithfully reproduces their curve. 


The calculation of u^., was based upon the wind data measured at the 

o 

18- and 30-m levels and the temperature data measured at the 18- and 60-m 
levels. Temperatures at the 30-m level were estimated by logarithmically 
interpolating between the 18- and 60-m levels. An estimate of the gradient 
Richardson number, equation (5), at the 23-m level (geometric mean height 
between the 18- and 30-m levels) was determined by assuming that the mean 
wind speed and temperature are logarithmically distributed between these 
levels. The gradient Richardson number estimated in this manner is given 
by 


iO 



-2 


, ( 12 ) 






g 



T(z,)-T(zi) g 


U(Z2) - U(Zl) 

T(z ) 

+ ^ 

1 z, C 

z In — ^ p 


z In 

g 

g Zl 


g Zl 


where T(z) is the mean temperature at height z, zj, and z^, denotes 18 

and 30 m, z = <\/z^~z^ 

S 

To calculate z was evaluated for each case by means of one 

of the three equations, (6) throu^ (8), corresponding to the appropriate 

Richardson number class. L* was then assumed to be invariant with hei^t, 

and iJr(18/L*) was estimated with equations (9) through ( 11) , Equation (3) 

was then evaluated at the 18-m level and solved to yield u, . The values of 

o 

used for this calculation are given in Figure 5. 

THE INERTIAL SUBRANGE AND REVISED VALUES OF THE 
SURFACE ROUGHNESS LENGTHS 

In the inertial subrange, the longitudinal spectrum is given by 


nS (n) 
u 


U' 


>:<o 





(13) 


where a is Kolmogorov's constant with a numerical value equal to 0. 146 
according to Record and Cramer [7] . The quantity <p^ is the dimensionless 

dissipation rate of turbulent kinetic energy per unit mass given by 


4 > 


e 


ki z e 
o 


(14) 


where e is the rate of dissipation of turbulent energy. Below 30 m, where 
the Monin and Obukhov similarity hypothesis for the wind profile is expected 
to be valid, <p^ is a function of Ri only. 
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Inferences concerning the dependence of (f>^ on Ri can be made with 

the aid of the eddy energy equation. For homogeneous terrain, this equation 
is given by 


dE 

dt 




A 

T 



(15) 


where — — is the eddy heat flux 0*w*, p is the mean density, E the 
O p 
p 

turbulent kinetic energy per unit mass, p* and w* denote the turbulent 
fluctuations of pressure and vertical velocity, and u , is the local friction 

velocity (u, = u, in the Monin layer) . Following Busch and Panof sky [8], 

h' 'I' O 

we write equation (15) in dimensionless form, so that 


kjZ dE 

u^l dt 
❖ o 




D 


(16) 


The terms in this expression are in one-to-one correspondence with ttiose of 
equation (15); however, the pressure term has been neglected. L is the 
Monin- Obukhov stability length where 


L 


u®, C p T 
*o p 

ki gH 


(17) 


Near ttie ground a majority of meteorological conditions are characterized, 
at least approximately, by horizontal homogeneity, steady mean wind with 
no change of wind direction with height, and steady heating from below. In 
most cases, it is reasonable to make these assumptions with regard to the 
I^C tower site. Thus, to a reasonable degree of approximation, we have 
dE/dt = 0, and equation (16) implies that in the Monin layer 




( 18 ) 
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Various authors have hypothesized schemes to balance the left-hand side of 
equation (18). Lumley and Panofslgr [6] suggest that the local mechanical 
energy production is balanced by the local viscous dissipation, so that 


= <P > ( 19 ) 

and thus, the buoyant energy production is balanced by the energy flux 
divergence term, making 





( 20 ) 


Busch and Panofsky [8] suggest the flux divergence term is negligible and the 
local viscous dissipation is balanced by both the local mechanical and buoyant 
energy productions, so that 




e 



( 21 ) 


We shall call equations (19) and (21) hypotheses I and II. 
According to Lumley and Panofsky, 


^ = (1-18 Ri)-^/'^ 


( 22 ) 


in the unstable Monin layer. This form of the dimensionless shear is con- 
sistent witii equation (6). Upon combining equations (13), (19), and (22), 
the spectrum in the Kolmogorov subrange for hypothesis I takes the form 


nS (n) 
u 

— 




(23) 
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Combination of equations (6), (13), (17), (21), and (22) yields 


nS (n) 
u 


u‘ 




- 2/3 
= a kj 


n 


K 


1 - Ri 


m 


(1 - 18 Ri)V^ 


2/3 


/ 


-2/3 


( 24 ) 


for hypothesis n, where and are the eddy viscosity and heat con- 
duction coefficients given by / 


K 

m 


du 

dz 


(25) 


and 


\ = - 


H 


C P~T~ 
p dz 


(26) 


Figure 6 illustrates ^nS^(n)/u?,^^J and |nS^(n)/u^. ^ J as functions 

-Ri for f = 1. 0 and Kj^/K^=1.3. As -Ri approaches infinity, 

|^nS^(n)/u^, converges to zero, and j^nS^(n)/u^^^J diverges. As -Ri 

approaches zero, jnS^(n)/u^ ^ j asymptotically approaches |^nS^ (n)/u? ^ j 

The difference between J and J is small (within the noise level of 


the data) for Ri > -1. 0, and thus it is difficult to reject one hypothesis in 
favor of the other for sufficiently small -Ri. 


At the 18-m level, longitudinal spectra were used to test hypotheses 
1 and II, and the cut-off value of f was approximately 2. 0 in most of the cases 
Actually, the Kolmogorov subrange occurs at much greater dimensionless 
wave numbers; however, the -5/3 power law behavior extends down to values 
of f on the order of unity for the longitudinal spectrum, while this is not true 
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Note: The variations of f nS (n)/u^, 1 for hypotheses 

L u 

I and II are indicated. The dashed curve corresponds 
to fnS (n)/u^, 1 with z = 0. 18 m and 

L f=1.0 ° 

= 0.23 m for a type I energy budget. 

Figure 6. The dimensionless logarithmic longitudinal spectrum at 
f = 1. 0 and z = 18 m as a function of -Ri. 

for the lateral spectrum. In addition, the Monin and Obukhov similarity 
hypothesis for the wind profile upon which the calculation of u ,_ ^ is based is, 

at most, valid below 30 meters. Accordingly, the 18-m longitudinal spectra 
are the only ones that could be used to test the validity of hypotheses I and II 
without introducing assumptions about how the eddy stress and heat flux vary 
with hei^t. 
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In Figure 6 we have plotted the experimental values of nS^(n)/u^ ^ 

for the longitudinal spectrum for f = 1. 0 as a function of -Ri. The experi- 
mental results scatter about the dashed line, but they appear to favor hypothesis 
I more than hypothesis H, especially for Ri :< - 3. 0. If we accept hypothesis 
I, then we must conclude that the scaling velocities, u^^ are too large and 

thus the roughness lengths are too large. If we reject hypothesis I, then we 
must accept a more complicated energy balance system. We shall invoke 
Ockham’s razor, accept hypothesis I, and correct the surface roughness 
lengths. 


We denote the surface roughness lengths in Figure 5 with z^^, and the 

correct ones resulting from the analysis of the longitudinal spectrum will be 
denoted with z^. The corresponding friction velocities will be denoted by 

u^^^^ and u^^^, respectively. The longitudinal spectrum in the inertial 

subrange, scaled in terms of u , , is given by 

oo 


nS^(n) 
❖ oo 


In -’i'(Ri) 

z 

oo 


In 


-'i^(Ri) 


a k^ (1 - 18 Ri) . (27) 


The data points in Figure 6 correspond to nS^(n)/u^ ^^, not nS^(n)/u? 

At neutral stability (Ri = 0) , we have ’t(O) = 0, and equation (27) reduces 
to 



for / = 1. 0. Denoting the right-hand side of equation (28) by ^d solving 

for z , we find 
o 


z 

o 


1-X X 
z z 


oo 


(29) 
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Upon extmpolating the data in Figure 6 to Ri = 0 at f = 1, we find 

nS (n)/u^ = 0.24, so that v = 1.059. The new values of z for this 

u oo o 

value of X z = 18 m are shown in Figure 7. Substitution of equation (29) 

into equation (27) yields 

. (30) 

o 

Zq (meters) 

1.0 
0.8 
0.6 
0.4 
0.2 
0 

0 90 180 270 360 

WIND DIRECTION (deg) 

Figure 7. Revised azimuthal distribution of the surface roughness 
length based upon energy budget considerations at z = 18 m. 

The dashed curve in Figure 6 represents nS^(n)/u^ as a function of Ri 
according to equation (30) for x = 059 and z^ = 0. 18 m (z^^ = 0. 23 m) , 

and it seems to fit the data reasonably well. The function nS (n)/ u^, 

for the range of variation of illustrated in Figure 7 departs from the dashed 



X ^ In — -'Jr(Ri) 
o 


In - 4r(Ri) 

z 


a (1 - 18 Ri) f 
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curve by only a few tenths of one percent. This means a spectral model of the 
longitudinal and lateral components of turbulence can be developed in terms of 
u^,^^^, and the final results can be corrected by applying a multiplicative factor 

which is a function of the Richardson number and a typical roughness length 
for the site. 


EXTRAPOLATION TO NEUTRAL WIND CONDITIONS (Ri = 0) 


The meteorological conditions of particidar engineering interest are 
those associated with mean wind speeds at the 18-m level which are greater 
than approximately 10 m sec“^. During these flow conditions, the boundary 
layer is well mixed so that vertical gradients of the mean flow entropy and thus 
potential temperature are small and the wind shears are large. Thus, the 
Richardson number vanishes or at least becomes very small. Accordingly, 
the neutral longitudinal and lateral spectra are of particular interest in the 
design and operation of space vehicles. The neutral spectra were determined 
by extrapolating the data to Ri = 0 by the procedure developed by Berman [ 9] . 
Scaled spectra nS(n)/u^, were plotted against Ri for various values of f, 

and smooth curves were drawn by eye. Of course, the data points scattered about 

this line. The values of nS(n)/u^, at Ri = 0 were then read off and 

*oo 

corrected by multiplying the results by to yield the neutral spectra 
nS(n)/u^, ^ for the various levels on the tower. The results of this graphical 

process are shown in Figures 8 and 9, where the positions of the maxima 
shift toward higher values of f as the height increases. This means that 
Monin coordinates j^nS (n)/u?, fj fail to collapse the spectra in the vertical 

so that F(f, Ri) is not a imiversal ftmction, and thus an added height 
dependence should be included in the analysis. Busch and Panofsky [8] have 
obtained similar results from analyses of tower data from Roimd Hill. The 
failure of the Monin coordinates to collapse the spectra in the vertical can be 
attributed to vertical variations in both the Reynolds stress and the length 
scale used to scale the wave number n/u(z). 

Above the Monin layer (z < 30 m) in the Ekman layer (z > 30 m) , the 
tangential Reynolds stress decreases with hei^t. In addition, the variances 
of the longitudinal and lateral components of turbulence are decreasing functions 
of z. Thus, if u^ is the correct scaling velocity, scaling the spectra with 

'I' 
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Figure 8. Dimensionless logarithmic longitudinal spectra for 
neutral wind conditions. 


the surface value of the friction velocity will cause the scaled spectra at the 
upper levels to fall below the 18-m spectra. 

By scaling the wave number with z, we have assumed that the integral 
scales of the longitudinal and lateral components of turbulence are proportional 
to z. One might suspect from the behavior of eddy coefficients [10] that, if 
the local integral scales have vertical variations, they should increase at a 
rate slower than z. In addition, we have no knowledge that the integral scales 
of the longitudinal and lateral spectra should have the same vertical variation. 
However, the analysis showed that Monin coordinates will collapse spectra 
with various tuitoulence intensities at any particular level in the vertical. 

To produce a vertical collapse of the data, it was assumed, for 
engineering purposes, that the spectra in Monin coordinates are shape-invariant 
in the vertical. This hypothesis seems to be reasonable and permits a 
practical approach to developing an engineering spectral model of turbulence. 
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Figure 9. Dimensionless logarithmic lateral spectra for neutral 

wind conditions. 

The Longitudinal Spectrum 

The vertical variation of the dimensionless wave number f associ- 

mu 

ated with the peak of the logarithmic spectrum scaled in Monin coordinates is 
given in Figure 10. A least-squares analysis of the data in this figure yielded 
the result 

where z is in meters. A plot of nS (n)/u^, versus f/f will shift the 

spectra at the various levels, so that all the peaks of the logarithmic longi- 
tudinal spectra are located at i/i = 1. Values of f from other tower 

mu mu 

sites are indicated in Figure 10. 
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Figure 10. Vertical distributions of the dimensionless frequencies 
f and f associated with the peaks of the logarithmic 

longitudinal and lateral spectra for neutral stability conditions. 

The average ratio /3^ of the shifted spectrum at level z and the 18-m 
spectrum, [\ is)] . is shown in Figure 11. A least- 
squares an^ysis of these data yielded the result 
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Figure 11. Vertical distributions of the collapsing factors /3 and 
for neutral stability conditions. 



where z is in meters. A plot of nS (n)/B u^, versus f/f will collapse 

u u =i'o mu 

the longitudinal spectra. The collapsed longitudinal data are plotted as a 

function of 0. 03 f/f in Figure 12. 

•"•'mu ® 

The function 


nS^(n) 

P u^, 

u *o 




mu 


1 + 1.5 


J- 




-1 5/3r. 


mu . 


(33) 
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0.03f/f_„ OR 0.1 f /f„„ 
mu mv 

Figure 12. Dimensionless logarithmic longitudinal and lateral spectra 

as functions of 0.03 f/f and 0.1 f/f for neutral stability 

mu mv 

conditions. 

was selected to represent the longitudinal spectrum, where and r^ are 

positive constants, determined by a least-squares analysis. For sufficiently 

small values of f, nS (n)/B u^, asymptotically behaves like f/f . This 

u ^u 'I'o mu 

behavior is correct for a one-dimensional spectrum. At large values of f, 

nS (n)/i3 u^, asymptotically behaves like (f/f \ consistent with 

u ^u \ muy 

the concept of the inertial subrange. The maximum value of equation (33) 

occurs at f = Various authors have suggested formulae like equation 

(33) to represent the longitudinal spectrum. However, most of the represen- 
tations have only one adjustable parameter available, while equation (33) has 
two: and r^. In this light equation (33) seems to be superior. 

controls the magnitude of the peak, r^ controls the peakedness, and f^^ 

determines the position of the peak of nS^(n)/u?,^^. Upon setting r^ = 5/3, 

we obtain the form of the longitudinal spectrum suggested by Panofslg?^ [6] to 
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represent the strong wind spectra of Davenport [11] . Von Karman's longitudi- 
nal spectrum [12] can be obtained by setting r^ = 2. A least-squares analysis 

of the longitudinal data in Figure 12 revealed that = 6. 198 and r^ = 0. 845. 


The Lateral Spectrum 

The lateral spectra can be collapsed with a procedure like the one 

used for the longitudinal spectra. However, to determine an analytical 
ejqtression for the lateral spectrum, special attention must be paid to the 
inertial subrange to guarantee that = 3/4 [13]. This requirement can 

be derived from the mass continuity equation for incompressible flow subject 
to the condition that the eddies are isotropic in the inertial subrange. The 
experimental values of and are given in Figures 10 and 11. These 

data show that and can be represented as power laws as for the 

longitudinal spectra. The function 


nS (n) 
v' ' 


f 


mv 


1 + 1.5 



(34) 


was used to represent the scaled spectra, where and r^ are positive 

constants. This function behaves like the one chosen for the longitudinal 
spectrum. 


For sufficiently large values of 
ratio between equations (33) and (34) 


f , the asymptotic behavior of the 
is given by 


S 


u 


S 

V 


c / f 

u "^u / mu 
C ^ \ / 

V *^v ' mv 


2/3 



(35) 
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In the inertial subrange we must have = 3/^> so that upon substituting 

this ratio into equation (35) , we obtain a relationship that can be used as a 
constraint in the determination of values of C and r and functions to 

V V 

represent B and f . The values C = 3, 954 and r =0.781, and the 
'^v mv V V 

functions 


f 


mv 



(36) 


and 



along with the longitudinal parameters, will satisfy the conditions of equation 
(35) and simultaneously give a good fit to the data (z is in meters). The 
collapsed lateral spectra and the functions given by equations ( 33)and (34) 
are shown in Figure 12. 


UNSTABLE SPECTRA 


To develop an engineering model for unstable conditions, the unstable 
spectra were averaged and then corrected by multiplication with 

r z 2 

In - (Ri) 

z 

o 

X“^ In — - 4^ (Ri) 

z 

o 


for z = 18 m, z^ = 0. 18 m, and Ri = -0. 3. The longitudinal and lateral 

spectra for the mean unstable conditions are shown in Figures 13 and 14. The 
unstable spectra were collapsed by using the procedures for the neutral bound- 
ary layer, and the functions given by equations (33) and (34) appear to be 
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equally valid for the unstable case. The functions f , f , B , B are 
^ ''mu "^mv '^u 

depicted in Figures 15 and 16, and the functions nS^(n)//5^u^^^ and 
nS^(n)/^^u^^^ are given in Figure 17. Table 1 summarizes the spectral 
properties of turbulence for unstable and neutral conditions. 



z (meters) 

Figure 15, Vertical distributions of the dimensionless frequencies 

f and / associated with the peaks of the logarithmic 
mu mv 

longitudinal and lateral spectra for unstable wind conditions. 

THE LONGITUDINAL AND LATERAL CORRELATION FUNCTIONS 

The normalized correlation function R(x) at space lag x is related 
to the spectrum through the Fourier integral 
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Figure 16. 


Vertical distributions of the collapsing factors 
and B for unstable wind conditions. 
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Figure 17. Dimensionless logarithmic longitudinal and lateral 

spectra as functions of 0. 04 f/f and 0. 033 f/f for 
^ •^mu •'mv 

unstable wind conditions. 





TABLE 1, SUMMARY OF THE SPECTRAL PROPERTIES OF 
TURBULENCE FOR NEUTRAL AND UNSTABLE CONDITIONS 


c ^ 

^ f 

nS(n) •'m 

m ^ 

^5/3r 


Neutral 

Unstable 

c 

6. 198 

2. 905 

U 



c 

3. 954 

4. 599 

V 



r 

0. 845 

1.235 

u 



r 

0.781 

1. 144 

V 



f 

‘^mu 

0.03 (z/18) 

0.04 (z/18)®-®^ 

f 

■’ mv 

0.1 (z/18)®-58 

0.033 (z/18)®*^2 


(z/18)"®*®3 

(z/i8)~®*^‘* 

^ 1 

(z/18) -“•3® 

ISI 

00 

1 

o 

k 


(t^R(x) =u J S(uk) cos (27Tkx) d« , (38) 

0 


where uS(uk) is the spectrum at wave number k (cycles m-^) and a is the 
standard deviation of the turbulence. The wave number is related to the 
frequency through Taylor’s hypothesis (k= n/u) . Substitution of equation 
(33) or (34) into equation (38) yields 


(t^R(^) 

^ O 




cos (27r4C)d^ 




5/3r 


(39) 
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where 


f 

f. 


\ 


> 


/ 


(40) 


The quantity ^ is the dimensionless space lag at height z, and the integral 
in equation (39) is a function of ^ only. Integration of equation (39) for 
neutral and unstable conditions with Simpson's rule yielded the results 
illustrated in Figures 18 and 19. 



Figure 18. Scaled correlation functions for the longitudinal and 
lateral components of turbulence as functions of the 
dimensionless space lag ^ neutral wind conditions. 

The dimensionless standard deviation can be obtained from 

equation (39) by setting ^ = 0 and then taking the positive square root of the 
resulting expression, so that 
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Figure 19. Scaled correlation functions for the longitudinal and lateral 
components of turbulence as functions of the dimensionless space 
lag I for unstable wind conditions. 
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jl/2^ 


'1' o 



d^ 


(l+ 1.5 ) 


5/3r 



(41) 


The right side of equation (41) is a pure number, values for which can be found 
in Table 2. 

The function 
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)3i^ 

=!'0 
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/3u- 
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(42) 
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TABLE 2. PROPERTIES OF THE CORRELATION FUNCTIONS 
FOR NEUTRAL AND UNSTABLE CONDITIONS 



Neutral 

Unstable 

a , 

U U 0 

2.227 

1.897 


1.677 

2. 302 


0.282 

0. 188 


0.332 

0.199 


was selected to represent the results of the numerical integrations for the 
neutral case. The parameter 6 is determined by least-square methods. For 
sufficiently small values of ^ , this function behaves like 


/3u2, 

^ ' o 




:2/3 


(1.5) 


(43) 


Now the theory of isotropic homogeneous turbulence predicts that, in the 
inertial subrange, 


R (x)-l 
V _ _4 

R (x)-l " 3 

u 


(44) 


and equation (43) predicts that 


R (x)-l 

V 

R (x)-l 
u 




2 

( 45 ) 


as The quantity within the brackets on the right side of equation 

(45) is equal to 4/3, according to equation (35), and 1. The 

apparent inconsistency between equations (44) and (45) results because 
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equation (44) is based upon the entire flow being isotropic, while equation (45) 
is based on spectra associated with turbulent flows which are only locally 
isotropic in wave number space for sufficiently large wave numbers. Thus, 
upon producing the Fourier integral, equation (39) , we obtained contributions 
to R(^) from both the isotropic and anisotropic portions of the turbulent flow. 
The quantities 


6=4. 758 
u 

6 = 3.399 

V 


and the function in equation (42 ) reproduce the results of the numerical 
integrations of equation (39) , for the neutral case, to within a few percent. 

In the unstable case, the function given by equation (42) does not 
reproduce the results of the numerical integrations at large values of | for all 
values of 6. To remedy this, an exponential factor was introduced into the 
expression and the function 




0 -^ 


1 + 6.815 


|3u^, 


(1.5) 


5/3r 


:2/s 



(46) 


was selected to represent the results of numerical integrations for the unstable 
case. The quantities X and y are determined by least-square methods. It 
appears that y = 0.9 can be used for both the longitudinal and lateral spectra 
and 


A =2.22 
u 

X = 2. 02 
v 

The function given by equation (46) yields a good fit near 4 = 0; however, it 
departs from the results of the numerical integrations by approximately 10 
percent at 4 = 1. Near the origin, equation (46) behaves like equation (43), 
consistent with hypothesis of the inertial subrange. 
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Let us now turn our attention to the longitudinal and lateral integral 
scales of turbulence which are defined by the expression 


CO 

L* = J R(x)dx . (47) 

o 


The integral scale, defined in terms of the dimensionless space lag, is given 
by 

L ' f °° 

^ = / R(^)d^ . (48) 


This integral is a pure number, and upon substituting the expressions given 
by equations (42) and (46) for R(4) and employing Simpson’s integration 
rule, we obtain the results in Table 2. 

Various investigators [4] have represented the correlation functions 
of atmospheric turbulence with the expression 


R(x) 



(49) 


The logarithmic spectrum associated with this expression is given by 


nS (n) 

=;<o 


27t/L* 

2 z 


(50) 


This function has a maximum value at 



z 


2tt 


0. 159 


(51) 
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This value of the dimensionless integral scale is significantly less (approxi- 
mately 45 percent) than the ones in Table 2 for the neutral case. However, for 
the unstable case, the integral scales given by equation (51) depart from the 
ones in Table 2 by only approximately 20 percent. 

In the neutral case, f and / are proportional to z and z®*®®, 

mu mv 

so that L* is a constant and L***. is proportional to z®*^^. The former is 
u V 

consistent with the results of Davenport's [11] analysis of high wind speed 

spectra. In unstable air, / and / are proportional to z®*®^ and z®*^®, 

mu mv 

and thus and L^‘ are proportional to z®*^® and z®*®®. 

THE DISSIPATION RATE OF TURBULENCE 


It is possible to estimate the dissipation rate of turbulence by 
examining the asymptote of the dimensionless logarithmic spectrum for large 
f values. According to equation (33), the longitudinal spectrum, for suffi- 
ciently large values of f , is asymptotically given by 


nS (n) 
u 

u^ 

O 


/3 


u 


f 


u 


(1.5) 


5/3r^j If. 


(52) 


mu 


Upon equating the right sides of equations (13) and (52), we obtain 


(1.5) 


/ C \ ®/2 

i 1 -JL ) gS/2 r 

5/2ru y Q! y ^u “^mu 


(53) 


According to the information given in Table 1 and equation (53) , 4> in 
neutral air is given by ^ 


0.055 


I 18 


(54) 
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while in unstable air, it is given by 



( 55 ) 


One should keep in mind that equation (55) was derived from a spectral model, 
which is an average of the unstable data, so that equation (55) is probably 
valid for conditions associated with values of Ri on the order of -0. 3. 

At z/L = 0 we have 0(0) = 1, so that, according to either hypothesis 
I or II [equations (19) or (21)], we must have 0^(0) = 1. However, equation 

(54) predicts that 0^(0) < 1 below z = 18 m for the neutral case and at 

z = 0 we have 0^(0) = 0. Actually, the layer of air in the domain 0 < z < 18 m 

is in the Monin layer, and we should have 0^ = 1 throughout this layer during 

neutral wind conditions. If we accept a 10 percent error in 0^ as a measure 

of the validity of the spectral representations given by equations (33) and (34) , 
then it follows from equation (54) that the model can be extrapolated down to 
the 3-m level in the neutral case. 

To understand the behavior of 0 in the unstable case, we write the 

€ 

eddy energy equation for the steady-state boundary layer in the form 





D 


= 0 


(56) 


Let us assume that the heat flux is height-invariant, so that L is a constant. 
At the surface of the earth, z/L vanishes, (u,/u, )^0(O) =1, and 

0j^(O) == 0 according to both hypotheses I and II; thus, 0^(0) = 1. As we 

proceed away from the earth into the unstable Monin layer, 0^^ is balanced 

by z/L, and 0^ is balanced by 0 according to hypothesis I (in the Monin 

layer u, =u, ). Because 0 is a decreasing function of -z/L, it follows 
'i< o 

that 0 is a decreasing function of z below 18 m. However, based upon an 
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analysis of a sample of Cape Kennedy data, Panofsky suggests that is 

unimportant in unstable air at 30 m and above. This conclusion may not be 
strictly true for the entire layer above 30 m, but it is reasonable to suppose 
that (p^ becomes negligibly small or vanishes somewhere above 30 m. In 

addition, for sufficiently large z, the dimensionless shear <p is small 
compared to -z/L and (u^,yu,^^)^ < 1, so that, according to equation (56), 

(p^ can be estimated as -z/L. Thus, at sufficiently large heights, cp^ is an 

increasing function of z. Therefore, <p^ must experience a minimum in the 

lower levels. Equation (55) implies that is an increasing function of z 

in the unstable boundary layer [Ri(18 m) ~ -0.3]; therefore, the minimum 
in must occur at the 18-m level if we accept hypothesis I and equation (55). 

The minimum in p^ actually occurs somewhere between the 18- and 30-m 

levels if the implications of hypothesis I are correct. If this minimum in p^ 

exists, it was not detected because of the wide spacing between the instrumen- 
tation levels on the tower. 

The level above 30 m at which p^ vanishes in the unstable model can 

be estimated with equation (56). Following Panofsky [14] , we assume that 
(u,/u^ )^p is small compared to -z/L; we set p =0 inequation (56) and 

then combine the resulting relationship with equation (55) to obtain 


= 0.63 


(t^) 


0.66 


(57) 


where z'*' is the level at which p^ vanishes. The 18-m level Richardson 

number for the mean unstable model is on the order of -0. 3. Therefore, 

according to equation (6) , 18/L* = -0. 19. Now, L = K L'/K, and 

m h 

K./K 1.3, so that L = -73 m. Substituting this value of L into 
n m 

equation (57) yields z*'* = 283 m. Above this level, p and p^ are small 

and p^ = -z/L. This should be compared with equation (55) , which predicts 

that p increases as z®*®®. 
e 
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The proposed energy budget scheme for the unstable boundary layer is 
summarized in Figure 20. In the lowest layer, the Monin layer, the energy 
budget is given by hypothesis I [equation (19)]. Above this layer, there is a 
transition region, between approximately the 18- and 30 -m levels, in which the 
energy budget transforms from a type I to a type III budget. A type III budget 
is one in which the mechanical and buoyant energy production terms and the 
energy flux divergence term all contribute to to varying degrees. From 

the base of the transition layer to a level on the order of 300 m, <p and (f>^ 

tend toward zero as z increases. Finally, above this level we have a region 
in which the buoyant energy production is balanced by the dissipation. For 
this scheme to work, the various functions in the energy equation must be 
functions of z/L and z scaled with a length scale other than L. The additional 
length scale should be a function of the external synoptic- and meso-scale 
conditions which force the turbulent flow in the boundary layer. 


VERTICAL VARIATION OF STRESS IN NEUTRAL AIR 


An estimate of the vertical variation of the stress in neutral air can be 
obtained with equation (54) and an inviscid estimate of e given by 


e 



(58) 


where is a representative dimension of the momentum carrying eddies. 

Jj 

In the Monin layer, 




(59) 


Above the Monin layer, increases slower than k^z and should approach 

a constant as z becomes large. In an analysis of the vertical distribution 
of wind in the neutral baroclinic boundary layer, Blackadar [10] points out 
that 


38 



r 


z (meters) 



Figure 20. Hypothesized scheme of the budget of turbulent energy 
in the unstable boundary layer for Ri(18 m) -0. 3. 


L 


B 


kjZ 

0.00027U 


(60) 


describes the few observed distributions of rather well, where U is the 

surface geostrophic wind and / is the Coriolis parameter. Upon combining 

c 

equations (54) , (58), and (60) , we find that, in the neutral boundary layer. 



(61) 
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This function vanishes at z = 0 and has a maximum value at 


z 

m 


18 


2.203 X 10~® 


U 

f 


( 62 ) 


where U and have units of m sec ^ and sec~^, respectively. At 30 deg 

north, / = 0. 728 • 10"^ sec“^; therefore, in the neutral boundary layer, 
c 

z /18 = 0. 03U. At z = 0 we should have (u /u, )^ = 1. However, the 
m V '■= o 

function given by equation (61) vanishes at z = 0 and increases with z when 

z < z . Thus, this function is valid only when z > z . The data which were 

used to construct the neutral spectral model were generally associated with 

conditions in which U < 15 m sec~^, making z < 8 m. Equation (61) pre- 

m 

diets reasonable estimates of the local stress for the netrual boundary layer 
above 18 m. The corresponding estimate for the unstable case with Blackadar's 
length scale fails to give reasonable vertical variations of the stress. 


CONCLUDING COMMENTS 

An engineering model of the longitudinal and lateral spectra of turbu- 
lence has been developed. The analytical expressions used to represent the 
spectra are asymptotically consistent with a sufficiently large wave number 
with the well-known properties of isotropic turbulence in the inertial subrange. 
Based on an analysis of the longitudinal spectrum in the inertial subrange at 
the 18-m level, it seems that in the Monin layer the mechanical and buoyant 
production rates of turbulent kinetic energy are balanced by the viscous 
dissipation and flux divergence terms, respectively. However, since the 
friction velocity was not measured directly, it is not possible to be sure that 
the correct hypothesis of the energy budget has really been selected. The 
friction velocities are derived for these cases from mean wind speed and 
temperature profiles by the use of assumptions whose validity can only be 
established if the friction velocity is known. Accordingly, the proposed energy 
budget should be considered tentative. 

The Fourier transform mates, or rather the correlation functions, 
associated with the spectral representations were obtained by numerical 
integration, and formulae were selected to represent the results of these 
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integrations. It is concluded from an analysis of the correlation functions 

that cr and a behave like and z~®*^^® in neutral air, while in 

u V 

unstable air they are proportional to and z~®*®^. In neutral air the 

longitudinal integral scale of turbulence is height-invariant, while the lateral 
integral scale behaves like In unstable air, L* and are propor- 

tional to z“®*^® and z~®*28, respectively. 

The vertical variation of the viscous dissipation was deduced from the 
spectral model. In the neutral boundary layer is proportional to z®*®®®, 

so that e is proportional to z”®*®^®. In the unstable boundary layer, (f>^ is 

proportional to z®*®®, so that e is proportional to z“®*®^. 

As one proceeds upward out of the Monin layer, the energy balance 
seems to become more complex, because the selective energy balance that 
occurs in the Monin layer fails in the Ekman layer and thus the mechanical 
and buoyant production terms together balance the viscous dissipation and flux 
divergence terms. At sufficiently great heights, on the order of 300 m, it 
seems that we again have a selective energy balance; however, we now have 
a situation in which the mechanical production and flux divergence terms are 
small and the buoyant production is balanced by the viscous dissipation. 

George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 

Marshall Space Fli^t Center, Alabama 35812, July 21, 1969 
933-50-02-00-62 
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